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(Dated: February 15, 2008) 

From an analysis of the flavor-tagged decay B® — > J/ip<ji we obtain the width difference between 
the B® light and heavy mass eigenstates, Ar s = Fl—Vh = 0.19±0.07(stat) Ip'Qi^yst) ps" 1 , and the 
CP-violating phase, <f) s = -0.571 9 ' 3o(stat) l^^syst). The allowed 90% C.L. intervals of Ar s and 
<f> a are 0.06 < Ar s < 0.30 ps _1 and —1.20 < 4> s < 0.06, respectively. The data sample corresponds 
to an integrated luminosity of 2.8 fb _1 accumulated with the DO detector at the Fermilab Tevatron 
collider. 



4 



PACS numbers: 13.25.Hw, ll.30.Er 



In the standard model (SM), the light (L) and heavy 
(H) mass eigenstates of the mixed B® system are ex- 
pected to have sizeable mass and decay width differ- 
ences: AA/ S = Mh ~ Ml and Ar s = T L - T H . The 
two mass eigenstates are expected to be almost pure 
CP eigenstates. The CP-violating mixing phase that 
appears in b — > ccs decays is predicted |l| to be <f> s = 
-20. = 2axg[-V tb V t */V cb V* s ] = -0.04 ±0.01, where V Z] 
are elements of the Cabibbo-Kobayashi-Maskawa quark- 
mixing matrix New phenomena may alter the phase 
to 4> a = -20 a + <£. 

In Ref . 0] , we presented an analysis of the decay chain 
B° s -> J/ij}(f>, J/i> -> fi+^r, (j> -> K+K~ based on 1.1 
fb _1 of data collected with the DO detector Q at the 
Fermilab Tevatron collider. In that analysis we mea- 
sured Ar s and the average lifetime of the B® system, 
t s = l/r s , where T s = (T H + T L )/2. The CP-violating 
phase 4>s was also extracted for the first time. The mea- 
surement correlated two solutions for <p s with two cor- 
responding solutions for Ar s . Improved precision was 
obtained by refittingthe results using additional exper- 
imental constraints [5(. Here we present new DO results 
of an analysis that includes information on the B® flavor 
at production time. Adding this information resolves the 
sign ambiguity on <fi s for a given Ar s and improves the 
precision of the measurement. The analysis is based on 
an increased data set, collected between October 2002 
and June 2007, and corresponding to an integrated lumi- 
nosity of 2.8 fb" 1 . 

We reconstruct the decay chain P° — ► J/tp<p, J/ip — * 
</> — > K + K~ from candidate (J/i(),(f>) pairs con- 
sistent with coming from a common vertex and having 
an invariant mass in the range 5.0 -5.8 GeV. The event 
selection follows that in Ref. [3| . The invariant mass dis- 
tribution of the 48047 candidates is shown in Fig. [T] The 
curves are projections of the maximum likelihood fit, de- 
scribed below. The fit assigns 1967±65 (stat) events to 
the P° decay. The flavor of the initial state of the P° 
candidate is determined by exploiting the properties of 
particles produced by the other b hadron ( "opposite-side 
tagging") and the properties of particles accompanying 
the P° meson ("same-side tagging"). The variables used 
to construct the opposite-side tagging are described in 
Ref. @ . The only difference to the description in Ref. Q 
is that the events that do not contain either the opposite 
lepton or the secondary vertex, and that were not used 
for the flavor tagging before, are now tagged with the 
event-charge variable defined in Ref. Q. 

Same-side tagging is based on the sign of an associ- 
ated charged kaon formed in the hadronization process. 
A B® (bs) meson is expected to be accompanied by a 
strange meson, e.g. K + (us) meson that can be used 
for flavor tagging. Such a configuration is formed when 



the initial b antiquark picks up an s quark from a virtual 
ss pair and the s antiquark becomes a constituent of an 
accompanying K + meson. Candidates for the associated 
kaon are all charged tracks with transverse momentum 
Pt > 500 MeV that are not used in the P° reconstruc- 
tion. We define the quantity AP = ^/(A<p) 2 + (A77) 2 , 
where A0 (A77) is the distance in the azimuthal angle 
(pseudorapidity) between the given track and the P s 
meson, and select the track with the minimum value 
of AP. The corresponding discriminating variable for 
the flavor tagging is defined as the product of the par- 
ticle charge and AP. Another discriminating variable 
is Qjet, the pr-weighted average of all track charges qi 
within the cone cos[Z(p, pb)] > 0.8 around the B meson: 

The discriminating variables of both the same-side and 
opposite-side tagging are combined using the likelihood- 
ratio method described in Ref. 0. A tag is defined 
for 99.7% of events. The performance of the com- 
bined tagging is taken from a Monte Carlo (MC) simu- 
lation of the P° — > J I ip4> process and is verified with the 
process for which we find the simulated 
tagging to be in agreement with data. The effective tag- 
ging power, as defined in Ref 0, is V = (4.68 ±0.54)%. 
It is a significant improvement over the performance of 
the opposite-side tagging alone, V = (2.48 ± 0.22)% 0. 
The purity of the flavor tag as a function of an over-all 
flavor discriminant is determined and parametrized, and 
the related probability P(B S ) of having a pure state P° 
at t = is used event-by-event in the fit described below. 



Data 
Total Fit 
Prompt Bkg 

non-Prompt Bkg 




5.5 5.6 5.7 5.8 
Mass (GeV) 

FIG. 1: The invariant mass distribution of the (J/xj),<j)) sys- 
tem for B° candidates. The curves are projections of the 
maximum likelihood fit (see text). 



We perform an unbinned maximum likelihood fit to the 
proper decay time, three decay angles characterizing the 
final state, and mass of the B® candidate. The likelihood 
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function C is given by: 



N 



(1) 



where N is the total number of events, and / s j g is the 
fraction of signal in the sample. The function !F^ ig de- 
scribes the distribution of the signal in mass, proper de- 
cay time, and the decay angles. For the signal mass dis- 
tribution, we use a Gaussian function with free mean 
and width. The proper decay time distribution of the 
L or H component of the signal is parametrized by an 
exponential convoluted with a Gaussian function. The 
width of the Gaussian is taken from the event-by-event 
estimate of the ct uncertainty <r(ct), scaled by an overall 
calibration factor determined from the fit to the prompt 
component of the background. J^ck ^ s ^ ne product of the 
background mass, proper decay time, and angular prob- 
ability density functions. Background is divided into two 
categories. "Prompt" background is due to directly pro- 
duced J/tp mesons accompanied by random tracks aris- 
ing from hadronization. This background is distinguished 
from "non-prompt" background, where the J/ip meson is 
a product of a B-hadron decay while the tracks forming 
the <p candidate emanate from a multibody decay of a B 
hadron or from hadronization. 

The decay amplitude of the B° s and B s mesons is 
decomposed into three independent components cor- 
responding to linear polarization states of the vector 
mesons J/tp and </>, which are either longitudinal (0) or 
transverse to their direction of motion, and parallel (||) or 
perpendicular (_L) to each other. The time evolution of 
the angular distribution of the decay products, expressed 
in terms of the magnitudes \Aq\, |, and and two 



relative strong phases Si 



of the amplitudes, is given in Ref. Q 



<5j_ and 82 



-S 



d 4 r 



dtd cos Odtpd cos tp 

2 cos 2 - sin 2 9 cos 2 ip) | A Q (t ) | 2 

+ sin 2 ip(l - sin 2 9 sin 2 tp^A^ 2 

+ sm 2 i;sm 2 6\Aj L (t)\ 2 

+ (l/v / 2) sin 2ip sin 2 9 sin 2tpRe(A* (t) A\\ (*)) 

+ (1/V2) sin 2tp sin 29 cos tplm(A* (t)A± (*)) 



- sin 2 ipsm29sm(plm(Al(t)A±(t)). 



(2) 



Polarization amplitudes for B® (upper sign) and B 
(lower sign) are given by the following equations: 



\A 0ill (t)\ 2 = IA,,n(0)| 5 
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sin <p s sin(AM s i) 



T_ =F e~ r * sin <p s sin(AM s i) 



ReK(t)A||(t)) = |^o(0)||A||(0)|oos(5 a 
sin(AM s t) , 
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lm{A\{t)A x {t)) = \A0)\\A x {Q)\ 



-Ti 



( ± sin 81 cos(AM s t) =p cos 81 sin(AM s f) cos (p s ) 



-(l/2)(e 



-Tut 



-r L t 



)sin0 s cos#i], 



where T± = (1/2) [(1 ± cos<p s )e~ TLt + (1 T cos <p s )er VHt ] 
For a given event, the decay rate is the sum of the B° s and 

B s rates weighted by P(B S ) and 1 — P(B S ), respectively, 
and by the detector acceptance. 

In the coordinate system of the J/ ip rest frame (where 
the (p meson moves in the x direction, the z axis is 
perpendicular to the decay plane of cp — > K + K~ . and 
p y (K + ) > 0), the transversity polar and azimuthal an- 
gles (9, ip) describe the direction of the fi + , and tp is the 
angle between p(K + ) and —p\J/tp) in the (p rest frame. 

We model the acceptance and resolution of the three 
angles by fits using polynomial functions, with param- 
eters determined using MC simulations. Events gener- 
ated uniformly in the three-angle space were processed 
through the standard GEANT-based 0] simulation of 
the DO detector, and reconstructed and selected as real 
data. Simulated events were reweightcd to match the 
kinematic distributions observed in the data. 

The proper decay time distribution shape of the back- 
ground is described as a sum of a prompt component, 
modeled as a Gaussian function centered at zero, and a 
non-prompt component. The non- prompt component is 
modeled as a superposition of one exponential for t < 
and two exponentials for t > 0, with free slopes and 
normalizations. The distributions of the backgrounds in 
mass, cos 9, tp, and cosi/' are parametrized by low-order 
polynomials. We also allow for a background term anal- 
ogous to the interference term of the A and A\\ waves, 
with one free coefficient. For each of the above back- 
ground functions we use two separate sets of parameters 
for the prompt and non- prompt components. 

The high degree of correlation between AM S , cp s , and 
the two CP-conserving strong phases 81 and 82 makes it 
difficult to obtain stable fits when all of them are allowed 
to vary freely. In the following, we fix AM S to 17.77±0.12 
ps -1 , as measured in Ref. [9(. The phases analogous to 
8i have been measured for the decay B® — * J/ipK* at the 
B factories. We allow the phases Si to vary around the 
the world-average values [10(] for the B^ — * J/ipK* decay, 
5\ = —0.46 and 8 2 = 2.92, under a Gaussian constraint. 
The width of the Gaussian, chosen to be 7r/5, allows for 
some degree of violation of the SU (3) symmetry relating 
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the two decay processes, while still effectively constrain- 
ing the signs of cos Si to agree with those of Ref. [Tpj ] . 
The mirror solution with cos S i < is disfavored on the- 
oretical 11] and experimental [l2| grounds. 



TABLE I: Summary of the likelihood fit results for three cases: 
free 4>a , <^s constrained to the SM value, and AV a constrained 
by the expected relation Arf M • | cos(<^ a )|. 





free (f> a 


^ = 0f « 


Arf 


T a (PS) 

Ar s (ps- 1 ) 

A ± (0) 

^(o)! 2 -]^!^)! 2 

5i 
S 2 

4>s 

AM a (ps^ 1 ) 


1.52±0.06 
0.19±0.07 
0.41±0.04 
0.34±0.05 
-0.52±0.42 
3.17±0.39 

-o w+ - 24 

u -° 1 -0.30 

= 17.77 


1.53T0.06 
0.14T0.07 
0.44T0.04 
0.35T0.04 
-0.48T0.45 
3.19T0.43 
= -0.04 
ee 17.77 


1.49T0.05 
0.083 ±0.018 
0.45T0.03 
0.33T0.04 
-0.47±0.42 
3.21T0.40 
-0.46 ±0.28 
ee 17.77 



Results of the fit are presented in Table |TJ The fit 
yields a likelihood maximum at 4> s = — 0.57^q'3q and 
Ar s = 0.19 ± 0.07 ps -1 , where the errors arc statistical 
only. As a result of the constraints on the phases Si, the 
second maximum, at S = 2.92^°;^, AT S = -0.19±0.07 
ps -1 , is disfavored by a likelihood ratio of 1:29. With- 
out the constraints on Si, (f> s shifts by only 0.02 for the 
AT S > solution. Confidence level contours in the S - 
AT s plane, and likelihood profiles as a function of <J) S and 
as a function of Ar s are shown in Fig. [2] Studies using 
pseudo-cxperiments with similar statistical sensitivity in- 
dicate no significant biases and show that the magnitudes 
of the statistical uncertainties are consistent with expec- 
tations. The mean value of the statistical uncertainty in 
<j) s from an ensemble generated with the same parame- 
ters as obtained in this analysis is 0.33. The test finds 
allowed ranges at the 90% C.L. of -1.20 < (j) s < 0.06 
and 0.06 < AT S < 0.30 ps" 1 . To quantify the level of 
agreement with the SM, we use pseudo-experiments with 
the "true" value of the parameter cj) s set to —0.04. We 
find the probability of 6.6% to obtain a fitted value of <f> s 
lower than —0.57. 

Setting (j) s = -2/3 a = -0.04, as predicted by the SM, 
we obtain Ar s = 0.14 ± 0.07 ps" 1 . This is consistent 
with the theoretical prediction of 0.088 ± 0.017 ps" 1 Q. 
The results for this fit arc shown in the second column in 



Tablcfl] The non-zero mixing phase is expected to reduce 
AT S by the factor of | cos(0 s )| compared to its SM value 
Arf M 0. In the third column of Table U we show results 
of a fit with AT S constrained by this expected behavior. 

The measurement uncertainties arc dominated by the 
limited statistics. Uncertainty in the acceptance as a 
function of the transversity angles is small, the largest 
effect is on |A (0)| 2 - |A||(0)| 2 . Effects of the imper- 
fect knowledge of the flavor-tagging purity are estimated 

^-0.4p 
'in 

Q. 

0.3 



0.1 
-0 
-0.1 
-0.2 



: (a) 




D0 , 2.8 fb" 1 
■ B°->J/\|/<|) 
AM S = 17.77 ps" 1 


_ / 1 




: — SM 

■ Ar = AT SM x|cos(<t> s )| 



1 1.5 
(radian) 




-3-2-10123 
(]> (radian) 



0.2 0.4 

Ar s (ps 1 ) 



FIG. 2: (a) Confidence-level contours in the AF S - 4> a plane. 
The curves correspond to expected C.L.= 68.3% (dashed) and 
90% (solid). The cross shows the best fit point and one- 
dimensional uncertainties. Also shown is the SM prediction, 
<p s = -2j3 s = -0.04, Ar 3 = 0.088 ± 0.017 ps" 1 Q. (b) 
Likelihood profile of <f> a , (c) likelihood profile of Ar s . 



by varying the flavor purity paramctrization within un- 
certainties. The "interference" term in the background 
model accounts for the collective effect of various physics 
processes. However, its presence may be partially due 
to detector acceptance effects. Therefore, we interpret 
the difference between fits with and without this term as 
a contribution to the systematic uncertainty associated 
with the background model. The main contributions to 
system atic uncertainties for the case of free <p s are listed 
in Table HJ 



In summary, from a fit to the time-dependent angu- 
lar distribution of the flavor-tagged decays B® — ► J/ilxj), 
we have measured the average lifetime of the B s ) 
system, t(B°) = 1.52 ± 0.05 ± 0.01 ps, the width dif- 
ference between the light and heavy B® eigenstates, 



Ar s = 0.19 ± 0.07(stat)+°;°i(syst) ps" 1 , and the CP- 
violating phase, 4>s = — 0-57i~Q3Q(stat)l"gQ2(syst). We 
also measure the magnitude of the decay amplitudes. In 
the fits, we set the oscillation frequency to AM S = 17.77 
ps -1 , as measured in Ref. ||, and we impose a Gaussian 
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TABLE II: Sources of systematic uncertainty in the results for the case of free <p a . 



Source 


T a (ps) 


AT 1 — 1 'I 
At s [ps ) 


A f(W 


\A (U)\ -\A\\[\})\ 




Acceptance 


±0.003 


±0.003 


±0.005 


±0.03 


±0.005 


Signal mass model 


-0.01 


±0.006 


-0.003 


-0.001 


-0.006 


Flavor purity estimate 


±0.001 


±0.001 


±0.001 


±0.001 


±0.01 


Background model 


±0.003 


±0.02 


-0.02 


-0.01 


±0.02 


AM S input 


±0.01 


±0.001 


±0.001 


±0.001 


±0.06, -0.01 


Total 


±0.01 


±0.02, -0.01 


±0.01,-0.02 


±0.03 


±0.07, -0.02 



constraint with a width of 7r/5 to the deviation of the 
strong phases from the values <5i = —0.46 and 5 2 = 2.92 
of Ref. {HI. The allowed 90% C.L. intervals of Ar s and of 
<j) s are 0.06 < AT S < 0.30 ps" 1 and -1.20 < (j) s < 0.06. 
The SM hypothesis for <fi s has a P-value of 6.6%. 

The results supersede our previous measurements @ 
that were based on the untagged decay — > J/ ty<$> and 
a smaller data sample. They are consistent with the re- 
cently submitted CDF results 
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